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1. Introduction 

Clarifying our understanding of B physics is an increasingly important undertaking in the hunt for 
Beyond the Standard Model (BSM) physics. The most recent round of experimental results from 
BaBar and th e LHC have imposed strict bounds on possible BSM contributions to rare B decays 
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |^ whilst continued experimental and theoretical progress has lead to more exacting 
tests of Cabibbo-Kobayashi-Maskawa (CKM) matrix unitarity SMIIinilS] 

Flavor-changing neutral current (FCNC) processes are forbidden at tree-level in the Standard 
Model of particle physics. Such processes can only proceed through loop contributions, which are 
sensitive to much higher energy scales than the b quark mass. FCNC decays therefore serve as 
a probe for BSM physics and can tightly constrain the nature and size of BSM interactions. In 
particular, precise measurements of the branching fractions of the rare — > decays offer 

two of the most promising avenues for the detection of BSM effects. These decays are strongly 
sensitive to the existence of BSM particles in a variety of BSM scenarios (see, for example, the 
discussion and references in^]). Semileptonic FCNC processes, such as B ^ K^*^ u'^ i i~ , rnay pro vide 
bounds on classes of BSM physics that are not ruled out by purely leptonic modes I51 14 | 15 | 16 | 

Taking a different tack, the B sector also provides the opportunity for detecting BSM physics 
through precision tests of the unitarity of the CKM matrix. This matrix, which is unitary in 
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the Standard Model, relates the mass and weak-interaction eigenstates of down-type quarks and 
incorporates all the flavor-changing and CP- violating couplings in the Standard Model. A wide array 
of channels probe the CKM matrix, such as the semileptonic decays B — >■ irtv and B — > D*lv, and 
neutral B mixing. Combining independent determinations of each CKM parameter over-constrains 
the elements of the CKM matrix and discrepancies or deviations from unitarity may indicate the 
breakdown of the CKM framework. Currently a number of tensions at the 2 — 3(7 level in unitarity 
fits hint at the exciting possibility of BSM physics I^U^. 

Whether the hunt for BSM physics involves FCNC processes or CKM unitarity tests, one thing 
is clear: lattice quantum chromodynamics (QCD) is an indispensable tool in the search. Precision 
tests of the Standard Model require both precise experimental and theoretical inputs. Quarks are 
confined to color-singlet states and the physics of the weak interactions underlying FCNC decays 
and CKM-probing processes must be teased out from the dominant nonperturbative dynamics of 
the strong force. Lattice QCD is one of the most important tools presently available for precisely 
determining nonperturbative effects, but until recently has, in many cases, been playing catch-up 
with experimental precision. Over the last five years, howev er, the era of precision la ttice QCD 
has dawned, with many results now at the few percent level [1111 8 | 19 | 20 | 21 | 22|23|24 | 2 5l the 
convergence of experimental and theoretical precision, and with growing experimental datasets and 
advancing lattice computations, tests of the Standard Model have become ever more exacting. 

What is less clear, perhaps, is how lattice perturbation theory fits into the march toward high 
precision tests of the Standard Model. What is lattice perturbation theory? And what exactly is 
the role of lattice perturbation theory in the hunt for BSM physics? This review will answer these 
questions. 

Viewed from the outside. Lattice QCD calculations are often seen as black boxes, simply non- 
perturbative number-crunchers that spit out ab initio predictions. Closer inspection reveals, not 
surprisingly, that this is not the whole story. I will attempt to lift the lid on heavy quark lattice 
QCD calculations — at least just a little — by examining the role of lattice perturbation theory 
in lattice B physics. I will frame this discussion in two ways: the first is the use of lattice pertur- 
bation theory in different formulations of heavy quarks on the lattice. Viewed through this lens, 
we see that lattice perturbation theory is a tool that is used, to varying degrees, by all the major 
lattice collaborations working on B physics. The second is the role of lattice perturbation theory 
in the calculation of specific physical processes, particularly key applications in the search for BSM 
physics. I focus on precision B physics; to this end I discuss not only processes currently in the 
headlights of BSM physics hunters, but also precision B spectroscopy, which has proved to be a 
trustworthy testing ground for lattice QCD and provides firm evidence that lattice QCD is indeed 
becoming more precise and reliable. For brevity I concentrate on the impact of lattice perturbation 
theory in lattice B physics, rather than discussing specific details of individual calculations. 

I start by motivating lattice perturbation theory in Sec.[5]and then briefly review recent methods 
for automating lattice perturbation theory in Sec. |31 In Sec. 2] I examine the ways in which lattice 
perturbation theory is used in various heavy quark formulations on the lattice. I briefly discuss 
relativistic formulations for heavy quarks on the lattice, but focus on the historically more prevalent 
effective field theory approach and use this discussion to limn the limits of lattice perturbation 
theory. In the second half of the review I highlight some applications of lattice perturbation theory 
to specific physical processes at the forefront of BSM searches and conclude with a short summary 
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in Sec. ID 

2. Lattice perturbation theory 

Lattice QCD is the preeminent approach for ah initio calculations of QCD processes and is conspicu- 
ously a tool for nonperturbative physics computations. In this view, the phrase "lattice perturbation 
theory" has the air of a contradiction — the perturbation theory of nonperturbative physics. But 
the phrase is slightly misleading: lattice perturbation theory is better viewed as "perturbation the- 
ory for lattice actions" . This move is not just cosmetic; the latter terminology is not only less 
outwardly troubling, but also generally more accurate and specific in scope. Unfortunately the new 
nomenclature is cumbersome; throughout this review I will continue to simply refer to lattice per- 
turbation theory (henceforth "LPT"), with the understanding that what is really at the heart of 
the discussion is "perturbation theory for lattice actions" . 

2.1. Motivating LPT 

A glance at the LPT literature reveals a range of applications within lattice QCD: calculating the 
renormalization parameters of bare lattice actions; matching regularization schemes and extracting 
continuum results from lattice data; and improving lattice actions. These applications, however, 
are inter-related and can be broadly categorized as accounting for the physics of energy scales 
excluded by the lattice cutoff. In other words, each of these applications falls under the banner of 
renormalization. 

Understanding these motivations for using LPT — that is, as a tool for renormalization — is also 
the key to understanding the justification for using lattice perturbation theory, laid out in^^SI. The 
lattice serves as an ultraviolet regulator discretizing spacetime and excluding all momenta greater 
than 7r/a (where a is the lattice spacing). As we would in any other regularization scheme, we must 
calculate the renormalization parameters of the regularized theory to correctly account for high 
energy effects. In this case the excluded scales are those above the lattice cutoff, which corresponds 
to approximately 5 GeV for current lattice spacings. At these energy scales, the coupling constant 
is relatively small, Q!s(vr/a) ~ 0.2, and perturbative approximations (and therefore LPT) are likely 
to be valid. 

LPT thus provides the connection between the low and high energy regimes of QCD. Moreover, 
this connection has been tested and validated in a wide range of QCD processes by comparing 
higher order perturbative calculations with nonperturbative computations in the weak coupling 

regime EZEHIMEIEaMl. 

There is, however, some subtlety to this story: we must choose our coupling constant carefully 
to ensure the perturbative series is well-behaved. Early LPT calculations were plagued by slow 
convergence and inconsistent results ^^ll These issues were the product of a poor choice of expansion 
parameter: the bare lattice coupling. We can greatly improve the convergence of ou r pe rturbative 
series by introducing an improved coupling constan t, ofte n defined in the "V-scheme" expressed 
at an appropriately chosen scale, the "BLM scale" I'^^ l '^^l . In some cases, LPT is simply insufficient 
and nonperturbative renormalization is required (see Section 14.21 for an example) . 

Before turning to look at some of the uses of LPT in different lattice formulations, I will examine 
in a little more detail both renormalization and improvement in the context of LPT. 
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Renormalization and matching calculations In general one has two choices for calculating 
renormalization parameters, matching parameters or improvement coefficients for lattice actions: 
nonperturbative tuning or perturbative calculation. There are a variety of nonperturbative tuning 
methods. A full review is beyond the scope of this report, but a partial list includes step-scaling 
methods, which may be applied to the Schrodinger functional off-shell Green functions (the 
"Rome-Southarnpton method" )'^^ or physical quantitieJSHl imposing Ward identities or chiral sym- 
metry relations 1-=^ and fixing to physical quantities, such as meson masses, or relations, such as 
relativistic dispersion relations 

Although these approaches differ considerably in practice, they share two particular advantages. 
All nonperturbative methods enable the replacement of perturbative truncation errors, which can 
be hard to quantify, with statistical and systematic errors, which can usually be determined reliably. 
In general, with intensified computational effort, these uncertainties can be systematically reduced. 
Secondly, a fully nonperturbative approach is the only truly ab initio method for QCD calculations. 

Nonperturbative methods share a number of disadvantages as well: step-scaling and iterative 
methods incur a greatly increased computational expense, whilst matching to physical quantities 
results in a loss of predictive power. 

LPT is an alternative to nonperturbative tuning. In this approach one calculates the renormal- 
ization constants or matching parameters perturbatively, carrying out the calculation in a spirit 
broadly similar to the standard perturbative calculations of continuum QCD. LPT is computa- 
tionally much cheaper than nonperturbative tuning and there is no loss of predictive power. There 
are disadvantages, too, of course. Most notably LPT introduces perturbative truncation errors. 
As lattice calculations become more precise, the concomitant perturbative error often become the 
largest source of uncertainty in the result, which will ultimately necessitate multi-loop calculations 
will eventually be necessary. This is no mean feat: LPT is usually more involved than the corre- 
sponding calculations in continuum QCD. The relatively recent advent of sophisticated automated 
perturbation routines for latt ice actions, however, means that multi-loop calculations for improved 
actions have been carried oufSSIllIllIlSl. Despite the complications, the computational cost of such 
calculations is invariably lower than that for nonperturbative tuning, ensuring LPT is an attractive 
alternative to nonperturbative renormalization. 



Improvement Precise lattice QCD computations generally require lattices with two properties: 
small lattice spacings to reduce discretization errors and large lattice volumes to remove finite size 
effects. Unfortunately lattice QCD calculations tend to be rather computationally expensive. And 
as the lattice spacing decreases and the lattice size increases, they grow ever more expensive. 

Thankfully for precision B physics there is an alternative to simply using ever finer lattice spac- 
ings. Lattice actions can instead be "improved" to eliminate sources of error. The most common 
approach to improvement is the Symanzik improvement programme '49t50!^ which scaling viola- 
tions are systematically removed by adding irrelevant operators to the lattice action, organized by 
mass dimension. This approach is frequently significantly cheaper than simulating at smaller lattice 
spacings. 

Determining the coefficients of the new terms in the improved action can usually be carried 
out quite simply at tree-level, but for precision studies, this is insufficient. By construction lattice 
actions exhibit the same long-range physics as continuum QCD, but their short-range behaviour is 
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distorted by the lattice cutoff. Radiative corrections to the new operators in the action renormalize 
the coefficients away from their tree-level values and this can lead to new sources of uncertainty in 
lattice predictions. By matching onshell quantities on the lattice and in the continuum, the resultant 
discrepancies can be removed order-by-order in perturbation theory. 

There is naturally a price to be paid for adding extra terms to the action (radiatively improved 
or otherwise). The extra operators considerably complicate the Feynman rules for the improved 
actions. Highly improved actions, such as the HISQ and NRQCD actions that I discuss in Sec. 21 
generate Feynman rules that cannot feasibly be manipulated by hand. Furthermore, the lack of 
Lorcntz symmetry complicates the evaluation of Feynman integrals, which are no longer amenable 
to Feynman parameter transformations or the other tricks of the continuum trade and must be 
handled numerically. To cope with these complications LPT has been automated. 

3. LPT today: Automation 

In recent years a number of automated lattice perturbation theory algorithms have been developed, 
many following the pioneering wo rk of Liis cher and Weisz"'-'^. An early variant of the Liischer-Weisz 
(LW) algorithm was deployed in | 42 | 52 | 53 | ^^-j- cm-j-ently the most widely used descendant of the 
LW algorithm is implemented in the HiPPy/HPsrc software packages l ^^ l ^^ l. More recently, the LW 
algorithm was adapted for gauge actions in the Schrddinger functional scheme A new software 
package, called pastor , has extended this work to include Wilson-t ype re lativistic quarks and 
HQET heavy quarks Initial calculations with pastor are underway l ^^ l ^^ l. Finally, following an 
altog ethe r different approach, a computer algebra system has been constructed and optimized for 
LPT EHl. Although the use of this software has so far been restricted to a calculation using the 
relativistic heavy quark action in the Columbia formulation ^ in principle the algorithm can be 
extended to any lattice quark action. 

The HiPPy /HPsrc softw are packages have now been used in a variety of perturbative calculations, 
for example inHSEHlTIHl 63 64 65^ ^^^^ extensively tested against previous results. Evaluating Feyn- 
man integrals is a two stage process with the HiPPy/HPsrc routines: one first generates the Feynman 
rules with HiPPy, a set of Python routines that encode the Feynman rules in "vertex files". These 
vertex files are then read by the HPsrc Fortran modules, which evaluate the corresponding integrals 
numerically via the Vegas algorithm or exact mode summation over a finite lattice. Any deriva- 
tives required in the calculation are irnplemented analytically using the derived taylor type, defined 
as part of the Fortran TayiUR package Both HiPPy and HPsrc implement parallel processing using 
MPI (Message Passing Interface). 

The Schrddinger functional scheme has been widely used over the last thirty years (see the 
Sec. 14.21 for references). The advent of LPT for the Schrddinger functional scheme is, however, 
much more recent. There are two reasons for this: the first is the additional complexity introduced 
by the Schrddinger functional formalism and the second is that LPT has very limited application 
to HQET computations, as I discuss in more detail in the next section. LPT is nevertheless useful, 
because perturbative calculations can serve as a starting point for more precise nonperturbative 
methods. These nonperturbative methods are computationally expensive and there are a large num- 
ber of renormalization parameters in the Schrddinger functional scheme, many of which have only a 
small relative effect on the final result. Perturbative investigations can guide where computational 
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resources should be focussed to most improve the precision of lattice calculations. 

There are three complications that must be handled by automated routines for the Schrodinger 
functional scheme. Firstly, translational invariance, which is assumed in the HiPPy/HPsrc algorithm, 
is broken in the time direction. Secondly, we must account for an induced abelian color background 
gauge fielcS. Thirdly, at a given order in the coupling constant, there are extra diagrams generated 
by the boundary conditions of the Schrodinger functional scheme that are not present in other lattice 
formulations. The pastor software package, written specifically to deal with these difficulties, is 
based on C++ routines that generate vertices and a Python frontend wrapper that allows the user to 
specify the lattice actions in symbolic form in C++. An xml input file specifies the desired observable 
in the Schrodinger functional scheme, which the Python frontend routine parses to generate all the 
relevant diagrams of o rder O n. These diagrams are then evaluated with a C++ program generated 
by the Python wrapper l ^^ l ^^ i 



The most recent automated LPT framework, presented based on a new C++ computer 

algebra library optimized for lattice perturbation theory. Lattice actions are defined by C++ classes 
that allow the user to specify the action in text form. These classes extract the vertices, then further 
C++ classes perform Wick contractions and convert the algebraic representation of the integrand into 
efficient C++ code that can be evaluated numerically. This framework also contains routines that can 
perform analytic differentiation with respect to the external momenta and und erta ke continuum 
calculations in dimensional regularization using the Passarino- Veltman reduction 



4. LPT today: Applications in heavy quark lattice formulations 

The computational expense of lattice QCD has been a distinct problem for b quarks, which are 
heavy and therefore "fall through" lattices that are too coarse to accurately resolve the Compton 
wavelength of the b quark. Putting this more precisely; it is not currently computationally feasible 
to use lattices large enough and fine enough to simultaneously satisfy rriTrL ^ 1 and mf,a ^ 1. Here 

is the mass of the lightest propagating particle in the theory; L is the side length of the lattice; 
and nib is the mass of the b quark. Heavy-light systems, which feature heavily in BSM searches in 
the B sector, represent the worst of all computational worlds; they are characterized by both the 
mass of the light quark and the b quark — scales separated by several orders of magnitude. 

Lattice QCD has traditionally dealt with this problem using the technology of effective field 
theories. Despite the considerable advances in computing infrastructure and simulation algorithms 
over the last t hirty years , even state-of-the-art lattice computations require extrapolations up to the 
b quark mass l ^3 | 46 | 47 |48| effective field theories continue to play an important part in lattice 
simulations. 

Historically, the application of effective field theoretic techniques to lattice QCD has elucidated 
the interplay of perturbative and nonperturbative physics and has consequently helped advance 
the understanding of conceptual issues such as the operator product expansion and the role of 
renormalons. I will examine effective theories for heavy quarks more closely to highlight some of 
the limitations of LPT, 

^In fact, the HiPPy/HPSRC packages have been extended to use background field gauge -S^, but this implementation 
is not sufficient for the Schrodinger functional scheme. 



November 28, 2012 2:4 WSPC/INSTRUCTION FILE mprreview 

Lattice Perturbation Theory and B Physics 7 



There are two widely used approaches to heavy quarks on the lattice: nonrelativistic QCD 
(NRQCD) and heavy quark effective theory (HQET). Both take advantage ol the large mass of 
the h quark relative to other scales in QCD physics and both have the same infinite mass, or 
stationary quark, limit. There are, however, important differences. Here I only briefly outline some 
of the important considerations for LPT for both lattice NRQCD and HQET. For a more complete 
pedagogical introduction to continuum NRQ CD, se e, for example, tSU. Reviews of HQET in the 
continuum appear in and in the textbooks whilst lattice HQET is reviewed in 

4.1. NRQCD 

Heavy quark bound states are typified by the small relative velocity of their constituent quarks. 
In heavy bb states, such as the r^b and T mesons, this relative velocity is approximately v'^ ~ 
0.1. Consequently, v induces three well-separated energy scales: the mass, 0{mb); the momentum, 
0{mbv)] and the kinetic energy 0{mi,v'^). NRQCD is also suitable for heavy-light systems, such as 
the B and Bg mesons, in which case the expansion parameter is Aqcd and the power counting 
is that of HQET, but for simplicity I restrict my discussion to heavy-heavy systems. 

To construct the NRQCD action we use the Foldy-Wouthuysen-Tani transformation to de- 
couple the quark and antiquark fields and undertake a nonrelativistic expansion. The result is an 
effective theory with nonrelativistic, low-energy degrees of freedom. At lowest order, the continuum 
NRQCD action is just the nonrelativistic Schrodinger action. 

Following the Symanzik improvement scheme, we can then systematically include interactions 
and improve the NRQCD action to a given order in v and Aqcd/^6 by adding irrelevant higher 
order operators. Finally, armed with an appropriately improved continuum NRQCD action, we then 
discretize it for lattice computations. We must choose the coefhcients of the higher order opera- 
tors to ensure the effective lattice theory replicates the correct behaviour for physical observables. 
As I discussed in the previous section, calculating the radiative corrections to the improvement 
coefficients is the role of LPT. 

One property of lattice NRQCD is a particular challenge: lattice NRQCD is a nonrenormalizable 
effective theory with no continuum limit. We cannot take the limit of vanishing bare lattice mass, 
arub — ^ 0. Even at leading order in the Foldy-Wouthuysen-Tani transformation, the 1/ {anib) correc- 
tions are present in the heavy quark propagator. Interaction and improvement terms in the NRQCD 
action are higher order in l/(amf,) and consequently there are an infinite number of divergences 
that cannot be absorbed into the renormalization parameters. 

From a perturbative standpoint, this means the series that define the radiative corrections to 
the higher order operators in the action become ill-defined at small masses. LPT — and indeed 
NRQCD — fails in the region arrib < 1. This has been explicitly demonstrated for the case of the 
chromomagnetic correction to NRQCD by studying the mass dependence for a range of masses 
around arub = 

lESl. We can safely use LPT to improve lattice NRQCD, provided we remain in the 
region for which amt, > 1. This ensures that nonperturbative computations using lattice NRQCD 
converge with those obtained from continuum QCD. 

Although we can not take the continuum limit of lattice NRQCD — it does not exist — we can 
nevertheless extract continuum results from lattice NRQCD calculations. By carrying out computa- 
tions at different values of the lattice spacing within the region of validity of NRQCD, i.e. at different 
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values of ami, for l/mt, < a < 1/Aqcd, we may parameterize the lattice spacing dependence of our 
results and quote uncertainties appropriately. With sufficient improvement, we can reduce the lat- 
tice spacing dependence to, for example 0{l/{ami,)'^) or higher, and obtain correspondingly precise 
results, often at the few percent level. Without LPT, such precise nonperturbative results for many 
B physics parameters would be a distant hope for the future. 



4.2. HQET 

Although HQET bears many similarities to NRQCD, the approach is both conceptually and prac- 
tically distinct. The first conceptual difference is that, unlike NRQCD, HQET is an effective theory 
for singularly heavy hadrons; mesons such as the T or r]i, are off limits. The second crucial distinc- 
tion is that HQET is believed to be renormalizable at leading order, i.e. in the infinite heavy quark 
mass, or "static", limit. The third contrast is that, even in the static limit, the renormalization and 
matchi ng parameters of lattice HQET must be calculated nonperturbatively for precise B physics 
resuhslJ^. LPT wiU not do. 

HQET is motivated by the intuition that singularly heavy hadrons are analogous to a hydrogen 
atom, with the b quark playing the part of the proton and the light quark the electron. Constructing 
HQET follows the same effective field theory approach to building NRQCD: we identify the rele- 
vant degrees of freedom, decouple the high and low energy modes via the Fold-Wouthuysen-Tani 
transformation and include all terms consistent with the desired symmetries at a given order in 
the expansion. We start by building an effective field theory that describes the light degrees of 
freedom interacting with a static color source, the heavy quark. Then, rather than including higher 
order operators directly in the HQET Lagrangian, we expand the non-leading terms in the path 
integral weight factor, exp(— S'hqet), so that the non-leading contributions appear as insertions in 
correlation functions. More explicitly, we write 



exp(-S'HQET) = cxp 



cxp 



a* 



(•^HQEt(2;) + •^HQEt(^)^ 

X 

E ^hq*^t(^)1 ( 1 - E 4'^et(^) ) , (1) 



where ^hqet is the leading order HQET Lagrangian and -C^ggrp contains higher order terms. This 
approach ensures that HQET remains renormalizable and the continuum limit of HQET correlation 
functions exists, provided all local operators at a given order in l/(am;,) are included. 

The parameters of the effective theory are fixed at tree-level by the FWT transformation, but 
beyond tree-level these must be found by matching to a l/mf, expansion of continuum QCD. Per- 
turbative matching is insufficient for precisi on B physics. Continuum HQET matching factors have 
been perturbatively studied in, for example, where it was demonstrated that even at three- loops 
the convergence of the perturbative series is very slow. The matching coefficients are scale depen- 
dent, and their magnitude can be reduced by lowering the scale, but this does not help: at the scale 
at which the perturbative series converges quickly the coupling constant itself is rather large. There 
seems to be no escaping the conclusion that for reliable results at the precision of a few percent, 
nonperturbative matching is necessary 
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Beyond leading order, matrix elements of higher dimensional operators mix and their coefficients 
must be fine tuned. These operators diverge with inverse powers of the lattice spacing and attempt- 
ing to remove these power divergences via perturbative matching introduces ambiguities associated 
with so-called "renormalons" . A full study of renormalons and their role in HQET would take us 
too far from our centr al topic; for in-depth reviews with applications to B physics and HQET, see 
for example, In essence renormalons arise because perturbative QCD is an asymptotic 

approximation to QCD that fails to capture nonperturbative behavior. A common example is the 
pole mass, which is a long distance quantity that can be rigorously defined order-by-order in per- 
turbation theory as the pole of the renormalized quark propagator, but is subject to ine scapable 
renormalon ambiguities of ©(AqcD/^if,) when nonperturbative effects are included I^HHSI, 

For our purposes, though, the central point is that renormalization parameters of HQET must 
be calculated nonperturbatively. This is a non-trivial task to implement for precision B physics, 
because one must ensure firstly that the HQET expansion is sufficiently accurate and secondly that 
the numerical precision is adequate. On top of this, the lattice simulations must be undertaken with 
physical volumes large enough to avoid significant finite volume corrections. Currently the ALPHA 
collaboration is carrying out a programme of precision B physics using HQET for the b quark 
| 24 | 25 | 84 | 85|86 | ^^-^^ implementing n onpert urbative renormalization using the Schrodinger functional 
scheme and a step-scaling method l ^^l^^l . 

All this does not mean, however, that LPT has no role to play in HQET computations. In fact, 
the ALPHA c ollabor ation has recently developed automated LPT routines for the Schrodinger func- 
tional scheme l ^^l^^l to serve as an exploratory guide to focus computational effort most efficiently 
in nonperturbative determinations. 



4.3. Relativistic heavy quark actions 

Relativistic heavy quark actions combine aspects of effective theories and ful ly relatiy istic treat- 
ments and hence serve as a suitable framework for both light and heavy quarks (89 | 90 | 91 | ^pj^g basic 
idea of the relativistic heavy quark formalism is to extend the Symanzik effective theory approach, 
on which NRQCD and HQET are based, to include interactions from both the small m^a and large 
mg/AqcD limits. As we take the limit of vanishing bare mass, m^a — ^ 0, the relativistic heavy 
quark action reduces to the O(a)-improved clover action Conversely, in the large bare mass 
limit, when ruq 3> Aqcd, the action breaks the time-space axis interchange symmetry and can be 
interpreted nonrelativistically, with a universal static limit. 

Both the Fermilab Lattice/MILC collaboration and th e RBC /UKQCD collaboration, currently 
employ relativistic heavy quarks for precision B physics l ^-^ l ^l l These collaborations' approaches 
to the relativistic heavy quark formalism differ primarily in the methods they use to tune their 
actions. 

The RBC/UKQCD method uses three physical conditions to tune their parameters, matching to 
the spin-averaged Bs mass and hyperfine splitting and ensuring the continuum dispersion relation 
holds for the Bg meson. This tuning process is nonperturbative, but there is nonetheless a role for 
LPT. Perturbative calculations of the three tuning parameters at one - loop h ave been carried out 
to provide a consistency check for the nonperturbative tuning process I I i Ultimately LPT will 
be used for parameters for which nonperturbative tuning is not available. 
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In the Fermilab approach, tw o of the parameters are fixed and the third tuned nonperturbatively, 
using the spin-averaged Bs mass^^. LPT played its part in the original construction of the Fermilab 
action and remains an important component for matching lattice matrix elements to their 
continuum counterparts. For decays such as the semileptonic B Triu and B Div decays, the 
Fermilab Lattice/MILC collaboration use a mixed approach to matching. The bulk of the matching 
is carried ou t nonperturbatively, leaving a factor close to unity that can be calculated perturbatively 
| 94 | 95 | 96 |97| pQj. other processes, such as neutral B mixing, however, the Fermilab Lattice/MILC 
collaboration uses a purely perturbative matching procedure 

Finally, it is worth noting that LPT played a c ritical role in the formulation of the relativistic 
heavy quark action used by the Tsukuba group ^^Hl this formulation, there is an extra parameter 
in the action that cannot be nonperturbatively determined ^^Sl. 



4.4. Fully relativistic b quarks: HISQ and twisted mass fermions 

Two lattice collaborations have recently begun to implemen t a progr amme of precision B physics 
using pur ely relativistic actions: the HPQCD collaboration | 46 | 47 | 48 | ^^^^ ^j^g ETM collaboration 
I23I99I100I 

These collaborations use different relativistic actions in their computations, but both presently 
require extrapolations up to the physica l b qu ark mass. The HPQCD collaboration uses the highly 
improved staggered quark (HISQ) actionH^, whilst the ETM collaboration employs a twisted mass 
variant of the Wilson action. The HISQ action exhibits exact chiral symmetry in the massless limit 
and therefore the heavy-light axial-vector and vector currents are absolutely normalized and do not 
require operator matching. The role of LPT in HISQ calculations is consequently somewhat reduced, 
but not completely absent. One loop calculations were required in the construction of the HISQ 
action. Furthremore, renormalization parameters for the tensor current and four quark operators 
are not absolutely normalized and therefore require matching to their continuum counterparts. 

The ETM collaboration use s no nperturbative renormalization for the quark bilinear operators 
used in B physics applications [iSSl Nevertheless, LPT is used to parameterize the discretization 
errors associa ted in the scale dependence of the nonperturbative renormalization factors in the RI- 
MOM scheme Subtracting the one- loop perturbati ve be haviour significantly improves the size 
of discretization errors in the renormalization constants ^IQ^. The part that LPT plays here may be 
restricted, but it has still been important in improving the precision of twisted mass computations. 

As above discussion makes clear, even as nonperturbative simulations move closer to the ideal of 
true ab initio calculations of QCD processes, LPT retains a vital role in extracting precise results. 
To clarify this role further, I now turn to the application of LPT in specific B physics processes. 



5. LPT today: Applications in B physics 

A handful of key B processes have become the loci of considerable experimental and theoretical 
attention. Effort has naturally focused on those processes that have, or are likely to soon have, 
small experimental and theoretical uncertainties. I will concentrate on only those channels that lie 
at the heart of the searches for BSM physics and CKM matrix unitarity violation and restrict my 
discussion to exclusive channels. Inclusive processes, whilst playing a vital role in BSM searches 
and contraining CKM matrix unitarity, fall outside the scope of this short review. 



November 28, 2012 2:4 WSPC/INSTRUCTION FILE mprreview 



Lattice Perturbation Theory and B Physics 11 

I summarize the processes in which LPT has been an important ingredient in Table [T] I break the 
table into three sections: in the top third I tabulate the key processes currently in the headlights of 
the flavor physics community. In the second third I present other processes for which LPT has been 
necessary for precise nonperturbative results. These processes have yet to achieve the experimental 
or theoretical importance of the processes listed in the first third of the table (usually because 
of larger experimental or theoretical uncertainties), but are nevertheless playing an increasingly 
important role in flavor physics in the heavy quark sector. In the final third, I present some more 
speculative approaches to uncovering new physics. Such processes are either poorly understood 
theoretically or are yet to be observed experimentally. 

Naturally this hierarchy is largely a difference of degree rather than difference in kind, but 
highlights the crucial processes at the center of current BSM searches. Within Table [1] I distin- 
guish between processes that determine CKM matrix elements, and thereby constrain CKM matrix 
unitarity, and rare decays that offer hope for detection of BSM loop effects. I represent this dis- 
tinction in the second column. In the third column I give the quantity determined by lattice QCD 
that characterizes the nonperturbative physics of the process. In column four I list the relevant 
renormalization parameters calculated using LPT. I use the shorthand Qq to represent heavy-light 
currents, QQ for heavy-heavy currents and QQqq to indicate four-fermion interactions. 



Table 1. Summary of B physics processes discussed in this article. For a full description of the table, 
see the accompanying text. 



Process 


B physics 


Lattice parameter 


Renormalization (LPT) 


Refs. 


B -i- ttIu 
B D*lv 
Bs — ^ Ai+A*" 
-Bq — B^ mixing 


rare decay 
VtdlVts 


form factors 
form factors 
decay constant 
5(7(3) breaking ratio 


Qq vector current 
QQ vector current 
Qq axial-vector current 
QQqq operator 


I62I97I115I123I 
I97I115I118I120I 

'97'115'123l 
,131il34TT35l 


B -5- Ti/ 

B Div 
B Ki+l- 
B ~¥ K'-y 


Vcb 
rare decay 
rare decay 


decay constant 
form factors 
form factors 
form factors 


Qq axial-vector current 
QQ vector current 
Qq vector current 
Qq tensor current 


I97I115I123I 
I97I115I118I120I 
I62I115I123I 
|62] 


B Xuiv 
Bs Ktv 

Be -> TV 


Vub 
Vub 

V,b 


form factors 
form factors 
decay constant 


Qq vector current 
Qq vector current 
QQ axial-vector current 


I62I97I115I123I 
I62I115I123I 
I97I115I118I 



5.1. Heavy -light current renormalization 

Matrix elements of quark currents are the starting point for any lattice computation of nonperturba- 
tive QCD contributions to weak interaction processes. Understanding quark current renormalization 
is consequently a critical component of precise lattice QCD predictions. 

The vector current mediates weak interactions between states of the same parity, such as the 
semileptonic decay B — ^ irtv, which is parameterized by 

{AV^\B) = U{q^) 



(2) 
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where /-|_ and /g are the "form factors". In the hmit — > 0, a good approximation for electron 
and muon neutrinos, the decay rate reduces to 

(3) 

This channel lies at the center of one of the main tensions in analyses of CKM matrix unitar- 
ity: currently, the values of Vub from inclusive and exclusive semileptonic decays exhibit a 3.3cr 
disagreement 

Other semileptonic B modes, such as B — > ptv, B — ?> rjlv and B uj£iy, represent rather 
more speculative possibilities for the extracti on o f Vub- The study of such decays is currently in 
its infancy for lattice QCD (see, for example, for a discussion of some of the challenges), but 
in principle any independent measurement offers the opportunity for greater understanding of the 
CKM mechanism. Similarly, the Bs ^(.v provides another exclusive determination of l yj) an d is 
curre ntly being explored on the lattice ^I^, in advance of experiments proposed at LHCb^I^ and 
Belled 

Purely leptonic decays of pseudoscalars, such as i? — >■ are directly sensitive to Vub and 
therefore offer further insight into this tension. Such decays require a spin-flip, however, an d the 
decay rate is helicity-suppressed by rnj/rn^. Only the B ^ tv mode has been observed l -^^^ i 
Consequently the experimental uncertainty of the branching ratio prevents determinations of Vub 
from reaching the precision achieved with the semileptonic modes. In addition, the leptonic decay 
B ^ TV is expected to be sensitive to the presence of a charged Higgs boson, which is an exciting 
opportunity for experimentalists, but means that determinations of Vub from this process may be 
difficult to interpret in terms of explicit unitarity violation. Nevertheless, experimental uncertainties 
will decrease and precise lattice computations will ensure that leptonic decays can ultimately deepen 
our understanding of the Vub discrepancy. 

Furthermore, the decay constants that parameterize leptonic decays are crucial ingredients in 
indirect searches for BSM physics. Observation of the rare — ?► decays is one of the key 

aims of the LHC. The high sensitivity of these processes to BSM physics ensures that they are 
amongst the star contenders for "likeliest channel for new physics results" . The branching fractions 
of these rare decays are proportional to the square of the decay constants, , defined via the 
matrix elements of the leptonic decays — ^ (.v as 

(0|Ao|B(,)) -/b,^„Mb(,), (4) 

where is the temporal component of the axial- vector current and Mb^,.^ the mass of the B^g) 
meson. Although the decay channel Bg Iv is not a physical process, it is nonetheless a process 
that can be computed within lattice QCD. In fact, the most precise values for the branching fraction 
of — >■ \i~ are obtained by taking a ratio with the 5° — B° mixing mass difference, AMg 
In this case, the dependence on js, drops out. Determinations using both methods have been 
remarkably consistent, providing a strong check of the reliability of lattice computations, because 
the lattice calculations for Bg decays and Bg — Bg mixing are rather different. 

Precise determinations of form factors and decay constants from lattice QCD are therefore 
of paramount importance in the hunt for BSM physics. And this is where LPT steps into the 
fray. The most up-to-date determinations of form factors and decay constants come from four 
lattice collaborations using four different formalisms. The ETM collaboration employ twisted mass 
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fermions with n/ = 2 species of sea quarks (with nj — 2 + 1 + 1 calculations underway) whilst 
early results from the AL PHA col laboration on n/ = 2 ensembles with_0(a)- improved Wilson 
fermions are now available EHH^. The Fermilab Lattice/MILC group IHMII and the HPQCD 
collaboration both use staggered fermion ensembles generated by the MILC collaboration, with 
different valence quark actions. The HPQCD collaboration has examined a wide range of flavored 
and unflavored meson decay constants. The current status of the decay constant "spectrum" from 
the HPQCD collaboration is illustrated in Fig. 15.11 and demonstrates the excellent agreement with 
experimental data. 
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All four collaborations use LPT, albeit to varying degrees. The HPQCD collaboration has made 
heavy use of LPT in recent computations, so I turn now to an example for the quark currents used 
by the HPQCD collaboration. 



5.1.1. A physics example: LPT in the HPQCD approach 

The HPQCD collaboration ha-s dete rmined both the decay constants Jb and /s^ | 47 | 19 | 11 5l 
the form factors /+ and /o I ^I^^ I H^ I. These determinations fall into two classes: the more tradi- 
tional calculations use nonrelativistic b quarks paired with relativistic staggered light quarks whilst 
recent computations have pioneered a completely relativistic treatment of the b quarks using the 
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HISQ action. This entirely relativistic approach has several advantages that simiiicantly improve 
the precision of the computation; the resulting value, /™^*^ = 0.225(4) GeVl^, is the most pre- 
cise currently available. Using relativistic b quarks avoids any need for effective actions and their 
associated systematic uncertainties, and the HISQ action has an absolutely normalized axial- vector 
current, removing the requirement for renormalization parameter calculations. 

This approach does not make LPT redundant, however. Simulations with relativistic h quarks are 
yet to reach the physical point and a determination of fs is even more computationally demanding. 
Lattice calculations with NRQCD remain the most competitive method for determining the ratio 
Systematic uncertainties associated with the perturbative matching or with neglecting 
higher order NRQCD terms are correlated between B and Bs simulations and therefore cancel 
in the ratio. The most precise result available for Jb takes advantage of these cancellations by 
combining a precise determination of the rat io / b / /b^ with the precise determination of using 
HISQ 6 quarks to obtain fs = 0.189(4) GeV^^^. This result would have not been possible without 
LPT. 



5.2. Heavy-heavy current renormalization 

Heavy-heavy currents are required as part of the extraction of the CKM matrix element Vcb via 
exclusive semileptonic decays. The decay B — ;> D*£iy is currently the most important exclusive 
process for attempts to reconcile the discrepancies between determinations of Vet from inclusive and 
exclusive channels, but Vcb can also be extracted from B — > Dliy. The BaBar collaboration recently 
reported an excess of 3. 4a ab ove the expected Standard Model ratio BR(_B — > _D'*Vi^)/BR(i? — > 
D^*^iv) in both channels ^^^Sl^ although a very recent lattice determination has reduced the tension 
to less than 2a 

Direct access to the CKM matrix element Vcb is possible through the leptonic decay Be — > tv^ 
which has decay rate 

T{Bc ^ TV) = ^Wbl'flml^^ f 1 - , (5) 

Stt "nig \ / 

where the decay constant is defined through 

(O|^o|i?(c)) = /B,.,MB,,,. (6) 

The experimental difficulties associated with reconstructing the r lepton from its decay products are 
considerable - and this decay is yet to be observed experimentally. This, however, offers a unique 
possibility for lattice theorists: the prediction of the leptonic decay rate from lattice computations. 
Successful predictions lend significantly more weight to our belief in the validity of lattice results 
than postdictions and are an important aspect of testing lattice QCD. 

Lattice determinations of the semileptonic B ^ D decays have so far been dominated by 
calculations from the Fermilab Lattice/MILC c ollab oration, which has studied various semileptonic 
B decays, including B ^> Dtv^, B D*lv^^Bxid the ratio oi~B D+i-V to B^ -> D+i'V 

Calculations of both the semileptonic B(^s) -D(s) and leptonic Be decays are underway by the 
HPQCD collaboration. 

Just as in the heavy-light sector, the role of LPT has been to determine matching coefficients for 
quark currents. The Fermilab Latticc/MILC collaboration's results for B ^ D processes required 
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the LPT heavy-heavy current matching calculations of The HPQCD collaboration is currently 
undertaking a determination of the matching coefficients required for the semileptonic B^^) — >■ 
and leptonic Be meson processes ^H^. 

5.3. B mixing 

The phenomenon of meson-antimeson oscillation is of particular interest to BSM physics hunters. 
Meson-antimeson mixing is a FCNC process and so, in the Standard Model, proceeds through one- 
loop box diagrams. This makes meson-antimeson mixing a particularly powerful contraint on both 
CKM unitarity and possible BSM physics. 

Meson-antimeson mixing occurs in both the B^ and B^ systems and is primarily mediated 
by top quarks. Determinations of the oscillation frequencies, IS.M and AM^ respectively, directly 
constrain the CKM unitarity triangle through the Standard Model relations 

AM, a \V:^Vt^\'flBB,, (7) 

where q may be either a d or s quark. Here the /s^ are the decay constants discussed in the previous 
section and the Bb^ are the renormalization group invariant bag parameters. The decay constants 
and bag parameters parameterize the nonperturbative matrix elements defined by 

(sJl (g7^(l - 75)6) {qr{l - Ts)^) = ^M|,/|,i?s„ (8) 

which are expressed in terms of the scale dependent bag parameters Bb^- The relation between 
BBg and Bb^ is known perturbatively to next-to-leading order ^ 8 | 129 | 

Although both AM and AMg directly constrain the apex of the CKM unitarity triangle, the 



SU(3) breaking ratio £, = jB^yBB^/fByBB^ provides a more stringent bound on CKM matrix 
unitarity because many systematic uncertainties cancel in this ratio. New physics may couple to the 
B'^ and systems differently and so separate constraints from both AM and AMg are required. 

Precise determinations of the matrix elements in eq. ([8]) require lattice QCD. And extracting 
continuum results from these lattice computations is once again a job for LPT. 

An exploratory study by the RBC/UKQCD collaboration found ^ — 1.13(12) usin g a single 
lattice spacing on n/ = 2 + 1 dynamical domain wall configurations with static b quark s LP T 
was used to improve the axial- vector current and match the results to the continuum | 131 | 132 | ^ 
more precise calculation by the HPQCD collaboration employed AsqTad light valence quarks and 
NRQCD b quarks on ensembles with nj = 2 + 1 AsqTad sea quarks lESl^ Using the perturbative 
matching results of ^'•^'^^ the HPQCD collaboration obtained = 1.258(33). These results are 
currently being updated using the HISQ action for the light valence quarks, which should reduce 
the uncertainty. Most recently, the Fermilab Lattice/MILC collaboration found ^ — 1.268(63) using 
the same ti/ = 2 + 1 AsqTad configurations, with AsqTad light valence quarks and b quarks 
implemented with the Fermilab action. Details of the perturbative calculation required for this 
result are forthcoming. 

5.4. b quark mass determinations 

Although not directly measurable, the mass of the b quark is nevertheless an important object of 
study, both as a free parameter of the Standard Model and as an input into searches for BSM 
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physics. 

The most precise nonperturbative resuh currently available for the h quark mass (expressed in 
the modified minimal subtraction, or MS^ scheme) is mjj^ {mj jq, n/ = 5) = 4.165(23) GeV, from a 
relativistic HISQ calculation by the HPQCD collaboration I^SI Two different determinations using 
Uf = 2 ensembles from the ALPHA and ETM collaborations are in agreement with the results 
from heavy HISQ simulations. Using nonperturbatively tuned HQET, the ALPHA collaboration 
obtained a preliminary result of 

^ATsi^Ais) = 4.22(10)(4)z GeV where the first error includes 
statistical and systematic uncertainties and the second arises from the quark mass renormalization. 
An alternative approach was taken by the ETM collaboration, using two different methods to find a 
(preliminary) average result of Wjgg(m;p;g) = 4.29(13)stat(4)sys GeV^^SI^ Here the first uncertainty 
is statistical and the second the total systematic uncertainty. 

5.4.1. A physics example: LPT in the HPQCD approach 

An example that illustrates the importance of LPT for precision B calculations is the extraction of 
the b quark mass from lattice NRQCD computations undertaken by the HPQCD collaboration. An 
early calculation obtained i^'M'sinT-'Ms) = 4.4(3) GeV^^, where the uncertainty was dominated by 
the o ne-lo op matching calculation needed to extract the continuum result. In recently completed 
work^l^, the one- loop calculation was extended to two loops using a mixed approach combining 
automated LPT and weak coupling simulations to extract the heavy quark energy shift. Preliminary 
indications suggest that this will reduce the uncertainties to the sub percent level: a real vindication 
of the power of LPT for precision lattice NRQCD. 

5.5. Heavy meson spectroscopy 

My focus has so far been the application of LPT to precision tests of the Standard Model in the 
hunt for BSM physics. These efforts require precise theoretical calculations, which, for nonpertur- 
bative QCD, often means lattice computations. But how can we be sure that the precision of our 
computations really is improving and that we are not missing sources of systematic uncertainty? 

To assuage these worries and to ensure that we really are improving the precision of our lattice 
calculations we need a proven testing ground for our results. In the B sector, this testing ground is 
heavy meson spectroscopy. 

The heavy meson spectrum is an excellent arena for testing and quantifying uncertainties asso- 
ciated with lattice calculations There are three reasons for this: firstly there are a number of 
gold-plated states that are experimentally and theoretically well defined; secondly there are mul- 
tiple states that can be used as parameter-free tests of the current status of lattice calculations 
and their associated uncertainties; and thirdly and most importantly, there remain a number of 
experimentally undiscovered states that allow lattice QCD predictions to be tested. Predictions 
have considerably more credibility that postdictions and provide the most stringent tests of lattice 
computations. 

An early vindication of the methods of precision lattice QCD occurred in the "prediction" of 
the Be meson mass. Although this mass was reconstructed with a precision of about 400 MeV 
by the CDF collaboration in 1998 the resolution was suffici ently poor that a later lattice 
calculation could claim a "prediction" of = 6304 ± 12(j"^* MeV^ i This result was confirmed 
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by observations at the Tevatron by the CDF collaboration L!^^ and the current world average value 
is me, = 6.277 ± 0.006 MevEl. The HPQCD recently updated its result to 6.280 ± 10 MeY^^, 
in excellent agreement with the experimental value. 

Heavy quark spectroscopy remains an important proving ground for lattice computations. The 
HPQCD collaboration has pioneered precision B spectroscopy using the HISQ action for light 
valence quarks and NRQCD b quarks. T his programme spans accurate postdictions for both heavy- 
light a nd h eavy-heavy systemsTI^IIlD! and p recis e predictions, such as the mass of the vector B* 
meson HS]^ the spectra for excited states qj. ^he D-wave spectra of the T The current 
status of the B spectrum from the HPQCD collaboration is illustrated in Fig. 15.51 
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The current status of precision B spectroscopy by the HPQCD collaboration. Figure reproduced courtesy 
Dowdall et al i*. 



The Fermilab Lattice/MILC collaboration has addressed heavy quark spectroscopy in 1 ^^ 1 143 1^ 
whilst the RBC/UKQCD collaboration has used relativistic heavy quarks to examine the low- lying 
heavy-heavy bound states The ETM and ALPHA collaborations are also currently undertaking 
calculations to address B spectroscopy. For heavy- heavy system s the spectrum and radiative decays 
of excited T states have recently been computed | 144 | 145 |146j -y^ork has also been carried out on 
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triply-heavy h bar yon spectroscopy, including not only the ground states I I I , but even the 
excited states 



5.5.1. A physics example: LPT in the HPQCD approach 

In the context of precision heavy meson spectroscopy, LPT has been primarily used to improve heavy 
quarks actions. The importance of this apparently small role is w ell illustrated by recent work on 
the hyperfine splitting of the T meson from lattice NRQCD In early determinations of this 
splitting, the HPQCD coUaboratiori used a tree-level value of the coefficient of the chromomagnetic 
interaction in the NRQCD action '^Sl. The results were not in good agreement with experimental 
data. A perturbative calculation of the one-loop correction to the chromo-magnetic interaction 
in NRQCD brought the lattice results in line with those from experiment This calculation 
demonstrates the effectiveness of the Symanzik improvement programme for lattice NRQCD and 
perfectly illustrates the importance of LPT to precision B physics. 

6. Conclusions 

The past decade has heralded a range of impressive successes for B physics on the lattice. Precise 
predictions and accurate postdictions, particularly in B spectroscopy, have convincingly demon- 
strated the effectiveness of contemporary lattice calculations. Improving determinations of decay 
constants and form factors have facilitated ever-tightening constraints on the unitarity of the CKM 
matrix and placed more stringent bounds on BSM physics effects in rare B decays. 

As more experimental data arrive, particularly from the LHC, further improvements and greater 
precision will be required from lattice computations. A close look at recent lattice QCD calculations 
reveals that lattice perturbation theory has been crucial in the development of precise theoretical 
results. A spectrum of applications of lattice perturbation theory can be found in contemporary 
lattice calculations, from the explicit — matching calculations used to extract continuum results — 
to the more implicit — in constructing, developing or improving lattice actions. Lattice perturbation 
theory has a home in the work of all the major lattice collaborations and will continue to help clarify 
our understanding of B physics. 
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